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ABSTRACT 

Q 

generalized force 


A high reduction ratio split torque gear train has been pro- 

q 

generalized coordinate 

pose< 

d as an alternative to a planetary configuration for the final 



stage of a helicopter transmission. A split torque design allows a 

r 

gear base radius 

high ratio of power-to-weight for the transmission. The design 
studied in this work includes a pivoting beam that acts to balance 
thrust loads produced by the helical gear meshes in each of two par- 

T 

total kinetic energy 

allel power paths. When the thrust loads are balanced, the torque is 
split evenly. A mathematical model was developed to study the 

t 

time 

dynamics of the system. The effects of time varying gear mesh stiff- 

V 

total potential energy 

ness, 

static transmission errors, and flexible bearing supports are 



included in the model. The model was demonstrated with a test 

X,Y,Z 

displacement in X,Y,Z direction 

case. 

Results show that although the gearbox has a symmetric con- 


figuration, the simulated dynamic behavior of the first and second 
compound gears are not the same. Also, results show that shaft 

• 

first time derivative 

location and mesh stiffness tuning are significant design parameters 
that influence the motions of the system. 

•• 

second time derivative 



a k 

angles defining locations of shafts (k = 1 

NOMENCLATURE 


helix angle of helical gears 

B 

energy dissipation function 

# 

angular displacement of gears 

C 

damping coefficient 

c 

damping ratio 

E 

mean value of displacement element function 

n 

fundamental mesh frequency 

e 

displacement element length 

(j 

input shaft speed 

K 

stiffness 


phase angle of time varying functions 

K 

mean value of stiffness 

Subscripts: 


L 

Lagrangian 

a 

axial 

M 

mass 

B 

bull gear 

m 

unit of length, m 

b 

bearing 

N 

unit of force, Newton 

bal 

balancing beam 


, 2, p, and B) 


1 



i input energy source 

j counter subscript for generalized coordinates 

a = 1,2, 3, ...18) 

o output load 

p pinion 

s shaft 

sB shaft between bull gear and output load 

sp shaft between pinion and input energy source 

sl,s2 shaft between spur and helical gear of 1st, 2nd com- 

pound gear 

xxbB bearing of bull gear in X direction 

xxbp bearing of pinion gear in X direction 

xxbl,xxb2 bearing of 1st, 2nd compound gear in X direction 

yybB bearing of bull gear in Y direction 

yybp bearing of pinion gear in Y direction 

yybl,yyb2 bearing of 1st, 2nd compound gear in Y direction 

zl,z2 axial direction of 1st, 2nd compound gear 

1,2 1st, 2nd compound gear 

lh,2h 1st, 2nd helical gear 

lmh,2mh 1st, 2nd helical mesh (pinion and compound gear) 
ls,2s 1st, 2nd spur mesh (compound and bull gear) 

INTRODUCTION 

Helicopter transmissions have evolved to a high degree of 
performance through the combined efforts of analysis, experiments, 
and application of field experience. F urther improvements in per- 
formance are desired, especially decreased weight, less noise, and 
improved reliability. Helicopter transmissions are typical of state- 
of-the-art systems in that computational tools are being used to de- 
velop and improve both new designs and existing components. 
Dynamic simulation of advanced transmissions can reduce develop- 
ment cost for new configurations compared to development by 
experimental methods which require expensive hardware. 

The most common final stage for a helicopter transmission is 
an epicyclic configuration which features an output shaft driven by 
several planets. This arrangement divides the transmitted torque 
among several planets. Epicyclic gear trains have been studied both 
analytically and experimentally (August et al., 1984; Choy et ah, 
1987; Choy et al., 1988; Zakrajsek, 1989). 

An alternative to an epicyclic stage which also transfers power 
to the output shaft through multiple paths is a split torque stage. 
One configuration for a split torque stage proposed for a helicopter 
application is shown in Fig. 1. White (1982) states that this new 
split torque gear train not only offers an overall weight reduction, 
but also promises the following advantages compared to the conven- 
tional epicyclic helicopter transmission: 

(1) High ratio of speed reduction at fined stage 

(2) Reduced number of speed reduction stages 



Figure 1.— Pictorial view of a split-torque helicopter transmission 
with two power paths. 


(3) Lower energy losses 

(4) Increased reliability of the separate drive paths 

(5) Fewer gears and bearings 

(6) Lower noise. 

The above advantages suggest the development of a new generation 
of helicopter transmissions of the split torque type shown in Fig. 1. 

A recent study on a split torque helicopter transmission is that 
of D. Hochmann et al. (1991). They analyzed the load distribution 
of spur and double helical gear pairs used in a split torque helicopter 
transmission. They concluded from their analysis that gear tooth 
profile modifications and staggered phasing of a double helical gear 
mesh can greatly reduce loaded static transmission error without 
seriously degrading load distribution. 

An aspect of split torque configurations that has not been 
rigorously studied is the overall dynamic behavior. A special char- 
acteristic of split torque drivetrains is the use of some method or 
mechanism to guarantee that the power is split evenly among the 
parallel paths. Several different methods have been proposed (Coy 
and Bill, 1988; Fisher, 1981; White, 1974, 1982, 1983, 1987, 1989). 
An appropriate analytical model is needed to study the system dy- 
namics of these designs. The purpose of this study is to provide a 
mathematical model to analyze the dynamic and vibration charac- 
teristics of split torque transmissions and to use the model to per- 
form some analytical studies. 

SPLIT TORQUE TRANSMISSION 

The split torque transmission considered in this study is shown 
pictorially in Fig. 1 and schematically in Fig. 2. The pinion transfers 
input torque to a set of two compound helical /spur gears via the 
helical meshes. The two spur pinions of these compound gears 
transfer the power into a final stage bull gear. The two helical 
meshes produce thrust loads that are reacted through a pivoted 
balancing beam. The function of the balancing beam is to couple 
the axial positions of the compound gears such that the thrust loads 
and therefore the torques of the two helical gears are equal. 

One of the main objectives of this work is to identify the key 
design parameters that influence the overall dynamic behavior of 
this system. Among the important parameters are the angles 
between the centerlines of the gears as identified in Fig. 2. These 
angles influence the stiffness properties of the system (appendix A). 
Also, for a given set of gears, these angles define the relative phasing 
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Figure 2. — Top view of the split-torque gear train with angles , 
a 2 , «p, and cxq depicting the shaft locations on the gear box. 


of the time-varying mesh stiffnesses. Note that although the gear 
train is geometrically symmetric with respect to the centerlines of 
the gearbox, the static bearing force reactions of the two identical 
compound gears are not equal. 

The mathematical model developed in this study may be used 
to understand and predict the effect of design parameters on the 
performance of the system. The important properties that can be 
simulated by the model can be categorized as: 

(1) Microgeometric parameters such as shaft hub eccentricities, 
tooth profile geometry, tooth spacing, lead, and profile errors. 

(2) Macrogeometric parameters such as gear mesh pitch, 
pressure angles, helix angle, and shaft and bearing locations. 

(3) Material properties and bulk dimensions of the components 
which defines the components’ inertia, stiffness, and 
dampening properties. 

(4) Type and geometry of the bearings that support the gear 
shafts. 


MATHEMATICAL MODEL 


The mass moment of inertia of the helical pinion, the mass 
moment of inertia of the balancing beam, and the mass of the 
balancing beam were considered negligible compared to the other 
inertias of the system. Under this assumption, the system is de- 
scribed by 18 equations of motion. The system is semidefinite and 
therefore has a rigid body mode. 

The equations of motion of this system were derived by the 
standard Lagrangian method as: 



dt 


£B 


dh 

dqj 


Qj (j = 1,2,3,... 18) 


a) 


where L = T - V, T = Total kinetic energy, V — Total potential 
energy, B = Energy dissipation function, q^ = Generalized coor- 
dinate, Qj = Generalized force associated with each generalized 
coordinate qj. All stiffness elements were assumed linear. The gear 
meshes stiffness elements were defined by time- varying periodic 
functions. 


A mathematical model has been developed to analyze the 
dynamic and vibration characteristics of split torque transmissions. 
The split torque transmission was modeled by a set of inertia, 
stiffness, damping, and displacement elements. The analytical 
model is shown in Fig. 3. Along with the inertia and stiffness 
elements shown in the figure, the model also includes a damping ele- 
ment parallel to each stiffness element, an input inertia, and an 
output inertia. Displacement elements are included at the gear 
mesh locations as illustrated in Fig. 4. The gear mesh stiffness and 
displacement elements together cam simulate the loaded static 
transmission error motions of the gear pairs. The displacement 
element is needed to simulate pitch errors and all other components 
of the static transmission error not attributable to the varying mesh 
stiffness. 

The mathematical model shown in Fig. 3 can be described by a 
total of 21 coordinates. The rotational positions of the gears, input 
and output inertias, and balancing beam require nine coordinates. 
The axial positions of the pinion, first compound gear, second com- 
pound gear, and balancing beam require four coordinates. The 
translations of the gear shafts require another eight coordinates. 


The kinetic energy of this system is: 

T = i £ Mjqf Qm 1,2,3,... 18) (2) 

The energy dissipation function of this system is: 

B = i £ C,q? (j = 1,2,3,... 18) (3) 

Here, the damping elements were modeled with proportional damp- 
ing expressed in terms of the damping ratio, £, stiffness coefficient, 
K, and inertia, M, as C = V 2£KM. 

The potential energy is classified into four groups of stored 
energy caused by: (l) distortion of the gear meshes; (2) shaft trans- 
lations at bearing supports; (3) twisting of gear shafts and hubs; 
and (4) axial distortions of the balancing beam components and 
support. 
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Figure 4. — Top view of the split-torque gear train, local and global reference frames. 


The first group is the stored energy due to distortion of the 
gear meshes. These terms were First expressed in a local coordinate 
system as shown in Fig. 4. Reference coordinate frames with super- 
script (*) are local. The potential energies specified in these 
coordinate frames were then expressed in a global reference frame 
via rotational coordinate transformations. These coordinate trans- 
formations introduce the shaft locations of the gear train as design/ 
analysis parameters. For example, the potential energy stored in 
the helical mesh of the input pinon and second compound gear was 
expressed in the local coordinate system as: 

V 2mh = 1 K 2mhM [ r p cos (^p - r 2h cos(/9) <t 2h + cos(£)Y* 

( 4 ) 

- cos(/J)Y 2 + sin(/J)Z p - sin(/J)Z 2 - e 2h (t) ] 2 
The local coordinate Y p was transformed by 


The potential energy stored in the other gear meshes was derived in 
a similar manner. 

The second group is the potential energy stored in the bearings 
due to shaft translations. To reduce the number of coordinates 
needed to describe the system, it was assumed that each gearshaft 
moves in pure translation and does not tilt. This assumption is 
reasonable when the length of the gearshaft is comparable to the 
face width of the gear. The bearing stiffness coefficients were 
defined in the global coordinate systems shown in Fig. 4. Under 
these conditions, the potential energy was expressed as: 


V b = \ K« bl X? + - 2 K yyb] Y 2 + i K Mb2 X 2 2 + i Y 2 


(V 


+ ^ K *xbp X P + \ Kyybp Vj + i K„ bB + i K yybB Y 2 


Y P = Y p cos ( a p) - V'^p) ^ 

By applying the relation 5 to Eq. (4), the potential energy can be 
expressed in global coordinates as: 

V 2mh = ^ K 2mh( C ) [ r p cos (^) 4 'p “ r 2h cos (^) 4, 2h + Coa(P) 

x cos(a p )Yp - cos(^)sin(a p )Xp - cos(/3)Y 2 ^ 


The third group is the potential energy stored in the twisted shafts, 
which was expressed as: 


V a = r K .1 f»!h - *J‘ 


X s2 i^2h 


( 8 ) 


+ 1 K sp i4-, - ¥ p j 2 + 2 K sBi * B - # 0 


+ sin (^)Z p - sin(0)Z 2 - e 2h (t)] 2 
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The fourth group is the potential energy due to distortion of the 
balancing beam support and the elements connecting the beam to 
the compound gear shafts. This was expressed as: 



v a = - K zl [Zi - L bal 0 bal - Z balJ 2 4- _ K 22 [Z 2 - L bal 0 bal 

2 2 (9) 

“ Z baiJ 2 + V 2 K bal Z bal + ^ Kz P Z p 

By applying the Lagrangian method, the equations of motion 
were derived. A typical equation of motion is presented in appen- 
dix A. Bearings were assumed to be pure stiffness elements with 
zero damping. The equations of motion are written such that the 
displacement element functions e^(t) appear on the right hand side 
of the equation as part of the generalized forcing functions Qj(t). 
The generalized forcing functions can be defined to simulate sources 
of excitation such as gear hub/shaft runouts and input/output 
torque fluctuations. However, in this study, only the static 
transmission error excitation as modeled by the time-varying 
displacement elements and gear mesh stiffnesses excited the system. 
The equations of motion were transformed to dimensionless forms 
using certain characteristic parameters inherent in the gear train 
system. Appendix B includes a glossary of the dimensionless 
parameters and how they were included in the mathematical model. 


ANALYSIS AND PARAMETRIC STUDY 

The mathematical model was demonstrated with a test case to 
simulate the motion of a split torque transmission. A gear mesh 
stiffness varies in time similar to a square wave function. The gear 
mesh stiffness elements were defined using a seven- term Fourier 
series. Seven terms were judged appropriate to approximate the 
square wave shape. A typical series used is illustrated in Fig. 5. 

The Fourier series functions used approximate the magnitude and 
wave shape of a typical gear mesh but are not precise represent- 
ations of a particular pair of gears. Displacement elements of the 
gear meshes were also defined using a seven- term Fourier series 
(Fig. 6). The displacement elements and stiffness elements together 
simulated typical loaded static transmission error motions. The 
fundamental period of each series was the corresponding gear mesh 
period. Dimensions and properties for geometry and inertias used in 
this study, given in appendix C, approximated those of the Split 
Torque Test Rig at the NASA Lewis Research Center. A damping 
ratio of { = 0.01 was assumed for all damping elements, and the 
bearings were assumed to be pure stiffness elements with zero 
damping. The method of R. Singh and T.C. Lim (1990) was used 
to calculate the bearing stiffness characteristics. 



Figure 5. — Periodic behavior of the stiffnesses of the helical 
meshes. 
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Figure 6. — Periodic behavior of displacement elements of the 
helical meshes. 

Parametric studies were completed to examine the effects of 
two design parameters on the performance of the gear train. The 
following parameters were studied: 

(1) The phase difference between the two helical gear meshes 
(stiffness tuning). There are two helical and two spur gear 
pairs in this gear system. Figure 5 illustrates a typical time 
variation of the stiffness properties. 

(2) The shaft locations relative to one another and relative to 
the frame of the gear box. Figure 2 shows that these relative 
locations can be characterized by the angles o^, a 2 , a p3 and 

The dimensionless equations of motion were integrated in time 
by a fifth-order/sixth-order Runge-Kutta method. 


RESULTS AND DISCUSSION 

The mathematical model developed in this work was used to 
simulate the vibratory behavior of the split torque transmission 
shown in Fig. 1. The source of excitation for the model was simu- 
lated static transmission errors for each gear pair. The stiffness 
element and displacement element at each gear mesh were defined 
appropriately to simulate typical static transmission errors for 
precision gears. All of the geometric and dynamic parameters were 
expressed in dimensionless forms. Figures 7 to 16 show some of the 
simulated motions of the system in generalized coordinates. These 
results were obtained for the case of an input pinion speed of 
u = 8778 rpm and an input pinion torque of 405 N-m. In this case 
tlfe pinion was positioned between the two compound gears such 
that the centerlines of the compound gears and the pinion lie on a 
straight line, in other words, at = 180°. 

The effects of relative shaft locations (a lf a 2 , a p , and a Q ) and 
of mesh stiffness tuning were studied. Here, mesh stiffness tuning is 
considered to be the adjustment of the phase difference between the 
periodic mesh stiffnesses of the first and second compound gears, 
^lmh “ ^2mh (Fig. 5). The results are contained in Tables I and II, 
respectively. Table 11(a) contains the maximum shaft orbit radii 
which are computed from X? + Z? for j = 1,2, p, and B. 
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Figure 7. — Displacement of the first compound gear shaft in 
X direction versus time. 
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Figure 8.— Displacement of the first compound gear shaft in 
Y direction versus time. 
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igure 9. — Displacement of the second compound gear shaft in 
X direction versus time. 
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Figure 10. — Displacement of the second compound gear shaft in 
Y direction versus time. 
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Figure 12. — Axial displacement of the second compound gear 
shaft in Z direction. 


Figure 14. — Dynamic transmission error of the pinion and second 
helical compound gear, 4i p - ij^h- 
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Dimensionless dynamic transmission error 



Figure 15. — Dynamic transmission error of the buil and first spur 
compound gear, 4* B - ^ 1s . 
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Figure 16. — Dynamic transmission error of the bull and second 
spur compound gear, i|/ B - ins- 


table I. - EFFECT OF SHAFT LOCATIONS ON GEARBOX 
VIBRATIONS 


(a) Lateral vibration of shafts 



Maximum dimensionless shaft orbit radius 

Shaft 

location, 

deg 

1st compound 
gear 

2nd compound 
gear 

Pinion 

Bull gear 

o, = 93.7 
a B = 50.0 

at p = 120.0 

0.4 

0.8 

B 

D 

Oj - 68.2 

a B as 57.3 
a p = 166.2 

0.5 

0.5 

B 

0.4 

- 61.1 
a B = 57.8 
a p = 180.0 

1.4 

n 


D 


(b) Axial vibration of shafts 



Maximum dimensionless vibration amplitude 

Shaft 

location, 

deg 

1st compound 
gear 

2nd compound 
gear 

Pinion 

a x = 93.7 
^ 50.0 
a p = 120.0 

2.5 

2.5 

0.3 

a t = 68.2 
= 57.3 
a p = 166.2 

2.0 

1.5 

0.02 

a, - 61.1 
a B = 57.8 
a p = 180.0 

2.0 

1.5 

0.02 


(c) Angular vibration of gear pairs 



Maximum dimensionless dynamic transmission error 

Shaft 

location, 

deg 

*p-*Lb 

-*ih 

mm 

a, = 93.7 
o B - 50.0 
a p = 120.0 

2.0 

1.0 

0.8 

a, - 68.2 
a B = 57.3 
a p s 166.2 

1.0 

1.5 

0.5 

a, ss 61.1 
a B = 57.8 
a p = 180.0 

1 

1.5 

0.5 













TABLE II. - EFFECT OF MESH STIFFNESS TUNING ON GEARBOX 
VIBRATION 


(a) Lateral vibration of shaft* 



Maximum dimensionless shaft orbit radius 

Phase 

difference, 

^lmh ” ^3mh' 

deg 

1st compound 
gear 

2nd compound 
gear 

Pinion 

Bull gear 

0.0 

0.4 

0.8 

1.2 

0.3 

3.6 

1.6 

1.1 

1.3 

0.3 

36.0 

16.2 

7.1 

1.7 

1.3 

90.0 

42.7 

17.1 

2.7 

3.0 

180.0 

53.65 

25.1 

3.6 

4.5 


(b) Axial vibration of shafts 



Maximum dimensionless vibration amplitude 

Phase 

difference, 

hnh ~~ 
deg 

1st compound 

gear 

2nd compound 
gear 

Pinion 

0.0 

2.5 

2.5 

0.3 

3.6 

2.5 

2.0 

0.2 

36.0 

2.5 

2.0 

0.2 

90.0 

2.0 

1.7 

0.2 

180.0 

1.8 

1.5 

0.2 


(e) Angular vibration of gear pair* 



Maximum dimensionless dynamic transmission error 

Phase 

difference, 


®P ” ^ab 


^ Ira fa " ^2rab, 

deg 




0.0 

2.0 

1.0 

0.8 

3.6 

2.0 

1.0 

7.0 

36.0 

4.0 

3.0 

60.0 

90.0 

5.0 

4.0 

160.0 

180.0 

10.0 

8.0 

200.0 


Table 1(b) shows the axial displacements of the first and second 
compound gears and of the pinion. Table 1(c) contains the maxi- 
mum dynamic transmission errors of the gear pairs. The numbers 
in these tables are given in dimensionless form. The physical 
displacements can be recovered by multiplying the dimensionless 
values by the characteristic parameters contained in appendix B. 

As shown in Table 1(a), the vibration amplitude of the first 
compound gear increases as the pinion is placed closer to the output 
bull gear; this closeness is measured by the angle a (Fig. 2). 
Installing the pinion closer to the output bull gear decreases the 
vibration amplitudes of the pinion shaft and has almost no effect on 
the output shaft vibrations. The data of Table 1(b) shows that 
bringing the pinion gear closer to the output gear reduces the axial 
movement of the compound gears and the pinion shaft. The data of 
Table 1(c) depicts that the dynamic transmission errors of the first 
compound gear mesh, ¥ — ^lh* an< ^ ^" irst out P u ^ s P ur niesh, 

decreases with moving the pinion toward the output gear 
while the dynamic transmission error of the pinion and the second 
compound gear, — ^ 2 h i ncreases - 

Tables 11(a) to (c) present similar results as those of 
Tables 1(a), to (c), except here, the mesh stiffness phase difference is 
varied to determine its influence on the response of the system. The 
data of Table 11(a) indicates that as the phase difference between 
the stiffnesses of the helical pinion and the two compound gears 


increases, the amplitude of lateral vibration increases for all shafts. 
However, an increase of the phase difference decreases the axial 
vibration for the pinion shaft and the two compound gears 
(Table 11(b)). The data of Table 11(c) indicates that an increase in 
the phase difference causes an increase in the dynamic transmission 
errors of the meshing gear pairs. These results show that the shaft 
locations and stiffness tuning are significant design parameters that 
influence the motions of the system. Optimization techniques 
should be applied with this analysis to determine the optimum shaft 
locations and mesh stiffness tuning parameters. 

The periodic, time-variable mesh stiffness categorizes the gear 
train as a parametrically excited dynamic system. The simulated 
motions of the system shown in Figs. 7 to 16 are consistent with 
typical responses of a parametrically excited system with high 
frequency characteristics restrained by a low-frequency envelope. 
Also, dynamic instabilities similar to those observed in systems 
expressed by Mathieu’s equations are expected under certain oper- 
ating conditions (Nayfeh, 1979). Unlike externally excited systems 
in which a small excitation produces large response only if the exci- 
tation frequency is close to one of the natural frequencies of the 
system, a small parametric excitation can potentially produce a 
large response when the excitation frequency is close to one- half of 
one of the natural frequencies of the system. In such cases, if the 
system is purely linear, theoretically the response amplitude 
increases to infinity. However, most physical dynamic systems 
possess some degree of nonlinearity which influences the system’s 
response as soon as the amplitude of motion grows beyond certain 
limits. Hence, the linear theory is only useful in determining the 
initial growth or decay of the system response. The dynamic system 
studied in this work has been modeled as a linear system with para- 
metric excitations. However, the physical system possesses some 
inherent nonlin ea r characteristics such as the gear mesh backlash 
and hardening bearing stiffnesses which influence the motion of the 
physical system. These nonlinear effects are not included in the 
theoretical model presented in this work. The large amplitudes of 
vibrations for $ 1# — as shown in Table 11(c) are thought to be 
the results of instabilities inherent in parametrically excited linear 
systems and not representative of the nonlinear physical system. 

As mentioned earlier, although the split-torque gear train is 
geometrically symmetric with respect to the centerline of the bull 
and pinion gears, the dynamic behavior of the first and second com- 
pound gears are not the same; this can be seen by comparing Figs. 7 
to 10, respectively. The difference between the vibrations of the two 
compound gears has been observed in an experiment conducted by 
the Sikorsky Helicopter Company as a part of the Advanced Rotor- 
craft Transmission Project (Kish, 1992). This difference may be the 
result of unequal static bearing reactions of the two compound 
gears. 

Figures 11 and 12 show the axial displacements of the two 
compound gears. The static transmission errors of the two power 
paths were considered slightly different in the dynamic simulation of 
the gear train; this is a realistic condition which will arise due to 
manufacturing limitations of making two exactly matched com- 
pound gear assemblies. The difference between the simulated static 
transmission errors may be the cause for the slight phase difference 
between the axial motion of the two compound gears, Z x and Z 2 , 
as depicted in Figs. 11 and 12. The phase shift between the axial 
motions of the two compound gears indicates that the balance beam 
has a rocking motion in the simulated dynamic operation of the gear 
train. 

Figures 13 to 16 are the dynamic transmission errors of pairs of 
the helical and spur gear meshes of the two power paths, respec- 
tively. The differences in the responses of the corresponding mesh 
pairs may be the result of the nonsymmetric bearing reactions of the 
power paths as well as the small differences in the simulated static 
transmission errors in each power path. 

The mathematical model presented in this work may be 
employed in an optimization process to determine the optimum 
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values of the geometric and material properties of the subcom- 
ponents of the gear train. 


SUMMARY OF RESULTS 

A mathematical model has been developed to analyze the 
dynamic and vibration characteristics of split-torque transmissions. 
The model consists of 18 equations of motion. The equations were 
developed using the Lagrangian method and solved using a Fifth- 
order/sixth-order Runge-Kutta method. The model was demon- 
strated with a test case. The source of excitation for the model was 
simulated static transmission errors for each gear pair. Parametric 
studies were also conducted to examine the effects of two design 
parameters, shaft location, and mesh stiffness tuning, on the per- 
formance of the gear train. The following specific results were 
obtained: 

1. The simulated motions of the sytem are consistent with 
typical responses of a parametrically excited sytem with high- 
frequency characteristics restrained by a low-frequency envelope. 

2. The large amplitudes of vibrations in some of the simulated 
motions are thought to be the results of instabilities inherent in 
parametrically excited linear systems and not represenative of the 
nonlinear physical system. 

3. Although the split-torque gear train is geometrically sym- 
metric, the simulated dynamic behavior of the first and second com- 
pound gears are not the same. 

4. The shaft location and mesh stiffness tuning are significant 
design parameters that influence the motions of the system. 
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APPENDIX A - TYPICAL EQUATION OF MOTION OF THE 
SPLIT-TORQUE TRANSMISSION 
The system has 18 equations of motion. 

The equation of lateral displacement of the first compound gear, Yl, 
is presented below: 

Mj Y, + {Cos 2 0?)K lmh + Cos 2 (a 1 )K lmj 
+ [r p Cos 2 (/3)K lmh ] 2 /[r p 2 Co 8 2 ( / 9)(K lmh + K 2mh ) + K, p ] 

+ KyyBl} Y 1 - {SinKJCosK) K lms^ ^1 

- {[( r p Cos 2 (£)K lmk )(r p Cc» 2 (/3)Sin(a p )K 2mJl )/[r p 2 Cos 2 (/?)(k lmh 

+ K 2mh ) + K. p ]} Xp 

- {Cos 2 (/J)K lmh 

- [(r p Cos 2 (^)K llnjl )(r p Co S 2 (^)(K lmh ^Cos(a p )K 2mh ))]/ 
[r p W(/J)(K lmh + K 2mh ) + K sp ;}Y p 

* {((r p Cos 2 (^)K lmh )(r p Co8 2 (^)K 2[nh i/[r p 2 Co3 2 (/3)(k lmh + K 2mh ) 

+ K sp } Y 2 - (CMa,) K lm ,} Y b 

+ {Cos^Sin^K,^ + [(r p Cos 2 (^)K ljnh )(r p Cos(/?)Sin(^)K lmh) i/ 
[r p 2 Co S 2 ( / 9)(k lmh + K 2mh ) + K sp ] } Z, 
-{[(r p Cos(^)Sin(^)K 2nlh )(r p Co S 2 (^)K 2mh ]/[r p 2 Cos 2 (^)(K Jmh 
+ IW) + K .pl} Z 2 - {Co9(yS)Sin(^)K lmll 
+ [(r p Co 8 2 ( / 9)K 1 ^)(r p Cos(/?)Sin( i S)(K lmh +K 2mh )]/ 

[r p 2 CosV)(K lmh + K 2mh ) + K„ p j}Z p 

- {[(K 9p )r p Co S 2 (^)K 2mh ]/[r p 2Cos 2 (^)(k lrah + ^2mh) 

+ {r lh Co S 2 (/9)K lmh +[(r p Cos 2 (^)K lmh )(r p r lh Co S 2 (^)K lmh )]/ 
[r p 2 Co S 2 (/9)(k lmh + K 2mh ) + K sp ]} * lh 
^{[(r p Cos 2 (^K lmh )(r p r 2h Cos 2 (^)K 2mh ))/[r p 2 Cos 2 (^)(K lmh 

+ K 2mh) + K S p]> *2h + { r uCos(a l )K llM } ¥ la 

- {r B Cos(a 1 )K lnM } ®|jj + C 2 2 Yj + C 21 X, 
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+ C 2i3 Xp t C 2i4 Y p + C 26 Y 2 + c 2>8 Y b 
+ C 2 g Zj + c 210 z 2 + c 211 Zp 

+ < -'2,12 + C 2j i3 *lh + *^2,14 ^2h + < -'2,15 *ls ~ ^2,17 

= - {[- K lmh Cos(£) + (r p 2 Co S 2 (^)K lmh 2 )/(r p 2 Cos 2 (^)(K lmh 

+ K 2mh) + VI e lh( l ) 

+ [(r p Cos 2 (/9) K lmh K 2mh )/(r p 2 Cos 2 (^)(K lmh + K 2mh ) + K, p )] 

X e 2h( l ) - l K lnu Cos (“l)] “ui 1 )} 

APPENDIX B - DIMENSIONLESS ANALYSIS 

The equations of motion are expressed in dimensionless forms, using 
certain characteristic parameters inherent in the gear train system. 
The characteristic parameters are: 

Length = E p , where E p is the mean value of the displacement 
element function of the pinion and helical compound gear 

Time = l/w p , where w p is the pinion speed, rad/sec 

Force = E p where K lmll is the mean value of the time 

varying mesh characteristic stiffness of the pinion-compound mesh 

Torque = (E K, mVl ) r , where r is the pinion’s base circle 
radius 

Angle = E p /r p 

Using the above characteristic parameters, the following 
dimensionless parameters may be defined: 

Time = wt 

Linear displacement = displacement /E p 

Angular displacement = (r p /E p ) angular displacement 

Inertia in linear displacement equation = M u p/K lmh 

Inertia in angular displacement equation — I w p /r p K 1 mh 

Damping_of linear displacement in linear equation of motion 
= ^ P C/K lmh 

Damping ofjinear displacement in angular equation of motion 
= «> p C/r p K lmh 

Damping of angular displacement in linear equation of motion 
= "p CAp Klmh 

Damping of ^angular displacement in angular equation of motion 

= "p CAp K lmh 


Stiffness of linear displacement in linear equation of motion 

= K(t)/R lmh 

Stiffness of linear displacement in angular equation of motion 
= K(t)/r p K lmh 

Stiffness of angular displacement in linear equation of motion 
= K(t)/r p K lmh 

Stiffness of angular displacement in angular equation of motion 
= K(t)/r p 2 K lmh 

Force = force/(E p K lmh ) 

Torque = torque/(r p E p K lmh ) 

The above dimensionless parameters were employed to carry out a 
computer simulation of the dynamics of this gear train. 

APPENDIX C - DIMENSIONS AND PROPERTIES OF SPLIT 
TORQUE TRANSMISSION 


DIMENSIONS AND PROPERTIES OF SALT TORQUE 
TRANSMISSION 


Gear design 

Helical 

pinion 

Helical gear 
of compound 
gear 

Spur pinion 
of compound 
gear 

Bull 

gear 

Number of teeth 

32 

124 

27 

176 

Pitch diameter (mm) 

51.1 

197.9 

68.6 

447.0 

Face width (mm) 

44.5 

38.1 

66.0 

59.9 

Normal pressure 

20.0 

20.0 

25.0 

25.0 

angle (deg) 
Helix angle (deg) 

6.0 

6.0 

.... 



Bearing design 

Input 

pinion 

roller 

bearing 

Input 

pinion 

ball 

bearing 

(Duplex) 

Compound 
gear roller 
bearing 

Bull 

gear bail 
bearing 

Bull gear 
roller 
bearing 

Inner race 
diameter (mm) 

50.0 

48.9 

87.4 

109.1 

113.0 

Outer race 
diameter (mm) 

69.1 

71.3 

66.5 

140.9 

133.9 

Number rolling 
elements 

13 

14 

15 

14 

23 

Rolling element 
diameter (mm) 

9.53 

11.13 

10.67 

15.88 

10.46 

Roller length 
(mm) 

13.20 


10.67 


10.41 

Contact angle 
(deg) 



29 


0 





Mass, 

kg 


Moment of 
inertia, 
kg-m 3 

Helical pinion 

2.0 

Input 

0.003 

Compound gear 

4.6 

Helical gear 

.010 

Bull gear 

29.5 

Spur pinion 

.003 



Bull gear 

.026 



Output 

.880 
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